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1. 


SPATIAL  SAMPLING  CRITERIA  FOR  NEARFIELD  MEASUREMENTS 
lyiADE  ON  CYLINDRICAL  SURFACES 

prA]is 

1 . This  report  establishes  the  spatial  sampling  criteria  for  nearfield 

measurements,  made  on  a cylindrical  surface  enclosing  the  sources,  that  enable 
the  farfield  of  those  sources  to  be  calculated  faithfully  from  the  nearfield 
data.  Examples  are  given  of  the  criterial  applied  to  both  omnidirectional  and 
directional  sources:  their  nearfield  and  calculated  farfield, 

CONCLUSIONS 

2.  ^The  sampling  criteria  for  nearfield  measurements  made  on  a cylindrical 
surface  have  been  established  in  this  report  for  both  axial  and  circumferential 
samples. 

3»  'c^.In  general  the  common  practice  of  using  half  wave  length  samples  circum- 
ferentially is  shown  to  be  adequate  for  omnidirectional  fields  (narrow  spatial 
bandwidths)  but  not  adequate  for  the  broad  circumferential  spatial  bandwidth 
fields  encountered  with  directional  sources.^  The  appropriate  sampling  criteria 
may  be  deduced  from  equations  (8)  and  (9)< 
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NOTATION 


^0  = ^0  Nearfield  co-ordinate  (see  Fig.  2), 


2L,  = 2Li  (p»y><^^ ) 


Farfield  co-ordinate  (see  Fig,  2), 


T(;^) 


The  acoustic  pressure  (amplitude  and  phase)  at  the  point  2^^  ■ 


T(....) 


The  Fourier  transform  of  the  acoustic  pressure  w.r.t,  <i>  or  z. 


T(....) 


w.r.t.  </>  and  z. 


dn^****' 


The  Foiu?ier  transform  of  the  normal  derivative  of  the  exact 
Green's  function. 


F.T. 

z < — * k ; k 


The  circumferential  spatial  wave  number. 


The  axial  spatial  wave  number. 


H (,...) 

u 

T (....) 

^ * * * * ^ 


Hankel  functions  of  the  first  kind  of  order  v. 


sampled  field. 


sampling  error  field. 


k = 27Tf/ 


the  temporal  wave  number. 


S (..,,) 


an  integer  ....  -1 ,0,+1 , ..,.  + 
The  Dirac  delta  function. 


i(7r(y+l)/2) 


~ H^(ka  sin  y) 


sampling  spacing  in  z, 
sample  spacing  in  <f>. 


INTRODUCTION 


4.  It  has  been  shown  by  several  authors  that  the  farfield  of  an  acoustic 
source  can  be  calculated  from  knowledge  of  the  acoustic  field  on  a specific  sur- 
face enclosing  the  source  (see  for  example  ref.  1»  2,  3).  This  calculation 
involves  the  evaluation  of  a surface  integral  of  the  pressure  field  of  the  source 
on  the  enclosing  surface. 

'■j.  In  the  case  of  a planar  surface  of  integration  this  surface  integral 
reduces  to  the  form  of  a Fourier  transform.  Nov;  the  acoustic  field  on  the  sur- 
face of  integration  can  only  be  determined  practically  at  discrete  points,  and 
the  required  spacing  of  these  samples  can  be  determined  by  applying  the  standard 
sampling  theory  to  the  field  and  its  Fourier  transform.  The  result  is  familiar 
that  the  samples  must  be  made  at  least  at  twice  the  highest  spatial  frequency 
(assuming  a spatially  band  limited  field  function  with  its  highest  spatial  fre- 
quency component  less  than  the  limit  of  propagating  waves). 

6.  In  this  report  the  equivalent  sampling  process  is  applied  to  a cylindrical 
surface  of  integration  and  the  corresponuing  sampling  criteria  established.  This 
IS  achieved  by  making  use  of  the  fact  that  the  surface  integral  relating  the  near 
and  farfields  can  be  considered  as  a convolution  integral  in  terms  of  the  cylin- 
drical polar  co-ordinate,  and  hence  use  may  be  made  of  the  convolution  integral 
for  Fourier  transforms  which  enables  the  sampling  problem  to  be  readily  analysed 
in  terms  of  the  Fourier  transform  of  the  cylindrical  field.  The  mathematical 
analysis  describes  both  axial  and  circumferential  sampling,  but  as  the  axial 
problem  is  essentially  the  same  as  for  a plane  surface  of  integration,  the  dis- 
cussion is  restricted  for  brevity  to  the  circumferential  sampling  criteria. 

7.  The  theory  applies  equally  to  nearfield  measurements  of  both  sources  and 
receivers. 

8.  Examples  are  given  of  the  sampling  criteria  applied  to  both  directional 
and  omnidirectional  fields. 


THEORY  (sec  Appendix) 


9.  The  farfield  of  an  acoustic  source  may  be  written  in  terms  of  its  near- 
field pressure  (ampl’itude  and  phase)  measured  on  a closed  surface  surrounding  the 
source  as,  (reference  1,3) 


) 


0Gf(r) 

«i!o)  — IS 


(1) 


where  T(X)  is  the  complex  pressure  at  the  point  X,  k,  is  a point  on  the  surface  S 
“ “ 9d(r) 

X^  is  a point  exterior  to  the  surface  S,  r = X^  ~ ""dn — normal 

derivative  of  the  exact  Green's  function  for  Dirichlet  boundary  conditions  on  the 
particular  surface  S,  with  respect  to  the  outward  normal. 

10.  Now  in  the  co-ordinates  of  figure  2 it  can  be  shown  (Appendix)  that  for  a 
cylindrical  surface  S with  Dirichlet  boundary  conditions  the  farfield  pressure 
IE  given  in  the  form  of  a convolution  in  0 of  the  neeirfield  pressure  with  the 
normal  derivative  of  the  exact  Green's  function,  i.e., 

T(a,s,0)  -2) 

- -r  ^ ~ ~ 

»jpREX;|QI^C  ’’4P^BLANK-NDT  FILMED 


T(p,y,0  ) 


= / 
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Relation  2 may  1)0  rewritten,  making  use  of  the  convolution  theorem  for  Fourier 

transformo  (reference  4)|  in  terms  of  the  Fourier  transforms  of 

3G(r) 

- ■■  With  respect  to  <f>  as, 

on 


1 1 ' - (i 

Tip,y,<P  ) = J ^ j — e 


do) . dz 
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vihere  is  the  Fourier  mate  of  and  represents  the  circumferential  spatial  wave 

number.  Now  by  using  the  particular  form  of  the  Green's  function  for  a cylindri- 
cal surface  with  Dirichlet  boundary  conditions  this  integral  may  be  re^^•ritten  as 
(see  Appendix). 

T(p,y,<^h  = ^~2 °v  y ^ (4) 

27t  a p ^ ^ 
v=«* 

t 

where  T Fourier  transform  of  the  nearfield  pressure  with  respect 

to  both  4>  and  z,  k is  the  Fourier  mate  of  z and  represents  the  spatial  wave 
z 

number  in  the  z axis. 

1 1 . The  errors  introauced  by  sampling  the  nearfield  pressure  can  now  be  found 

» 

by  substituting  for  T (k^,w^)  by  the  wave  number  spectrum  of  the  sampled  near- 
field which  then  allows  relation  4 to  be  expanded  into  a term  v^Ihlch  gives  the 
correct  field  and  an  error  term,  (see  Appendix) 

T’(.PtY><p‘' ) = T(p,y,<^^ ) + Tg(p,y,</)^ ) (5) 

where  T(p,y,0”*)  is  the  correct  field  given  by  relation  4,  and  Tg(p,y,<#)'*  ) is  the 
error  term  given  by. 


v»,y.4)  = i I y y 

^ 271  ap<t>z  ^ ^ 

O®  I p=-o*  q=-e 


^ T(k^  - 27rp/_,  - 277q/-) 


(q/o) 


Thus  the  farfield  calculated  from  the  sampled  nearfield  has  been  separated  into 
the  correct  field,  and  an  error  term  which  contains  the  true  wave  number  spectrum 
of  the  nearfield  repeated  about  points  2ttp/-  and  2vq/-  in  and  respectively. 

Now  by  a suitable  choice  of  the  sample  spacings  z and  ^ these  repeated  spectra 
may  be  removed  from  the  ranges  of  and  k^  that  contribute  to  the  farfield  thus 

allowing  the  correct  field  to  be  calculated  from  the  sampled  nearfield.  This  is 
discussed  further  in  the  next  section* 


r 


7. 


DISCUSSION 


L 


12.  It  is  shown  above  that  the  farfield  of  a source  can  be  described  by  the 
inverse  Fourier  transforx  of,  the  product  of  the  transform  of  the  neai’field 
pressure  on  the  cylindrical  surface  with  the  transform  of  the  normal  derivative 
of  the  exact  Green's  function  for  the  surface.  The  normal  derivative  of  the 
exact  Green's  function  for  a cylindrical  surface  with  Dirichlet  boundary  condi- 
tions is  shovm  in  Fig.  3a, b together  with  its  Fourier  transform  with  respect  to 
(j>  calculated  from  I'elation  (A4).  The  modulus  of  the  Fourier  transformed  deriva- 
tive of  the  Green's  function  decays  rapidly  beyond  some  value  of  as  can  be 

seen  in  Fig.  3b.  The  region  of  for  which  the  modulus  o1'  the  transform  is 

effectively  non  zero  will  depend  on  the  accuracy  required  for  the  particular 
problem  and  will  also  depend  on  the  value  of  ka  sin  y.  This  is  shown  in  Fig.  3c, 
For  most  practical  purposes  values  of  the  normalised  modulus, 


, 3G(w  ) 

I ■ I 

I 0n 


3G(o) 

3n 


> 10 


-10 


are  considered  non  zero  and  the  corresponding  region  of  u,  can  be  seen  from  Fig  3c 
to  be  well  approximated  by  the  relation 


$ M ka  sin  y + C 


where  M = 1 and  C = 1 6 in  the  range,  10  < ka  sin  y < 50.  In  the  limit  of  large 
ka  sin  y the  gradient  of  the  line  tends  to  1 . Should  greater  accuracy  be  required 
then  the  approximate  value  of  M and  C can  be  determined  from  Fig  3c.  Now  the 
Fourier  transform  of  the  nearfield  pressure  on  the  cylindrical  surface  has  been 
shown  due  to  the  .'•ampl  ing,  to  be  repeated  about  points  in  spaced  ^tt/-  apart  (see 

Reference  5,  Fig  and  Appendix),  The  region  of  this  transformed  field  that  con- 
tributes to  the  farfield  is  that  part  which  is  enclosed  by  the  effectively  non 
zero  region  of  the  modulus  of  the  Fourier  transformed  derivative  of  the  Green's 
function, 

I ! 

I I 

Thus  the  sampling  interval  (p  must  be  small  enough  to  ensure  that  the  repeated 
spectra  do  not  overlap  into  this  region.  This  may  be  r-aimncr' 3<  t a;;  ’ p 


2t7 


<t> 

where  ui 


> w . M ka  sin  y + C 

<pmax 


(?a) 


pmax 


is  the  highest  circumferential  spatial  wave  number  of  the  nearfield 


pressure.  Rearranging  equation  7a  gives  the  required  angular  sample  spacing  <p  as, 

^ < 2vj 


(M  ka  sin  y + C + 


and  the  circiunferential  sample  spacing  as, 
X < K I 


{ sin  y sin  y ]{ 


(7b) 


(8) 
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Consider  now  two  special  cases 

Broad  Circumferential  Spatial  bandwidth 

1 3.  The  nearfield  pressure  in  this  example  is  taken  as, 

T(z,0)  = cos^  «xp  (i  2va{^  - cos  <^)/a.) 

which  is  assumed  for  the  purpose  of  analysis  to  be  typical  of  the  field  of  a 
directional  source.  This  function  and  its  Fourier  transform  are  shown  in 
Fig  3a, b for  ka  sin  y = 25.8, 

14.  The  maxim’um  spatial  frequency  in  this  example  can  be  taken  approximately 
as  50.  Hence  from  reiation  7 the  angular  sample  spacing  is  given  as  4>  < 0,06 
radians  or  in  terms  of  the  circumferential  samples  x < a/(4  sin  y). 

15.  The  effect  of  various  sampling  inteirvals  on  the  predicted  farfield  is 
shown  in  Fig  oa,b  together  with  the  spectra  of  the  nearfield  for  ka  sin  y = 25.8. 
It  can  be  seen  that  a/4  samples  will  indeed  allow  the  farfield  to  be  calculated 
with  an  error  term  less  than  about  95  dlB  below  the  pattern  maximum;  somewhat 
mor-e  accurate  than  normally  required^  However  with  nearfield  samples  every  A/2 
spacing  the  error  term  in  this  example  is  only  25  dB  below  the  pattern  maximum. 

So  although  a sampling  criterion  of  A/2  spacing  is  in  this  case  shown  to  be 
inadequate  a criterion  of  a/4  spacing  is  certainly  sufficient. 

Narrow  circumferential  spatial  bandwidth 

16.  The  nearfield  pi’essure  in  this  example  is  taken  as 

T(z,d.)  = 1 

which  is  an  ideal  non-directional  circumferential  field.  The  maximum  circum- 
ferential spatial  wave  number  is  obvious,  max  = 0,  Hence  from  relation  76 

and  taking  ka  sin  y = 25.8,  M =1.4  and  C = l6  we  have  that  the  circumferential 

sampling  criterion  is  x < a/(2  sin  y). 

1 /.  The  efl'ect  of  various  sampling  intervals  on  the  calculated  farfield  is 
shown  in  Fig  7a, b together  with  tlie  nearfield  spectra.  It  can  be  seen  that  the 
nearfield  is  faithfully  calculated  with  A/2  samples;  with  0,8a  samples  the  error 
is  about  0.5  dB  and  with  1 .bA  samples  about  8 dB. 

Circumferential  and  Axial  spatial  samp'J  ing  criteria 

18.  It  can  be  seen  above  that  the  error  term,  relation  b,  will  be  negligible 
if  the  f ollowin,-'  inequalities  are  satisfied,  (see  Fig  4).  For  axial  samples. 


2tt  , 

— > k + X 


zmax 


• e.  z < 


Yo<  ♦ 


zmax 


(9) 


.vher-c  k = — — , and  is  the  limit  of  real 
c ’ 

axial  spatial  wave  number  present.  This 
critei'ion  when  k„  max  = k. 


angles  for  y,  and  k is  the  highest 
^ ' ’ zmax 

leads  to  the  familiar  a/2  sampling 
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For  circiiml'erential  sanpling, 

^ ^ (M  ka  sin  y + C + (7b) 

or 

^ ^ Jsin  y [M  + (C  + y]l  (®) 

wherf-  w.  is  the  highest  c ircumierential  spatial  wave  number  and  M =1,4  and 
(pmax 

C = 16  for  most  practical  purposes.  Kquations  (8)  and  (9)  are  shovm  evaluated 
in  figs  ba,b. 

CONCLUSIONS 


iy.  The  sampling  criteria  for  nearfield  measurements  made  on  a cylindrical 
surface  have  been  established  in  this  report  i’or  both  axial  and  circumfei'ential 
samples. 

20.  In  general  the  common  practice  of  using  half  wave  length  samples  circum- 
ferentially is  shov/n  to  be  dequate  for  omnidirectional  fields  (narrow  spatial 
tandwidths)  but  not  adequate  for  the  broad  circumferential  spatial  band-width 
fields  encountered  with  directional  sources.  The  appropriate  sampling  criteria 
may  be  deduced  from  equations  (8)  and  (9). 
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AITKNDIX  1 : MATHEM.ATICAL  DETAILS 

1 . The  I’ield  at  a point  ^ exterior  to  a closed  surface  3 containing  the 
sources  may  be  written  for  Dirichlet  boundary  conditions  on  the  surface  S. 

• . 9G(r) 

9G(r) 

where  T(X„)  is  the  field  at  a point  on  S and  — r is  the  normal  derivative 

'—0  9n 

of  the  Green's  function  for  Dirichlet  boundary  conditions  on  S. 

2,  For  an  ini'inite  cylindrical  surface  S y/e  may  write  (see  Fig  2j 


9G(r) 


i!c(p  - z cos  y)  +« 


^ 2 
2n  ap 


1 . . i(n  (v  + 1 ;/2; 
H (.ka  sin 


or  for  simplicity 


9G(r) 


(z.y)  2_  °‘u  i<P  - <f>  )') 


where  (z,y)  = 


ik(p  - z cos  y)  i(m  (v  + 1 )/2) 

5 and  a (y)  = — — jr ^ r— 

„ 2 y H (ka  sin  y j 

27T  ap  y ' ' ^ 


Equation  A1  may  now  be  rewi-itten. 


• • 1 

-r,/  ' ■■■ml  j.\  9G-  (a,z,p,y,<f>  - d>  ) j . , 

T(p,y,^  ; 3 , T (a,z,0)  — ' d 0 d z 


v/hich  has  the  form  of  a convolution  in  0 and  hence  we  may  Fourier  transform 

T (a,z,0)  and  (a>ZiP,y,0  0 ) respect  to  0 and  use  the  convolution 

theorem  for  Fourier  transforms  to  give, 

1 • • . 9G  (zjW  ) ia),0 

T(p,y,0  ) = J T (z,m^)  ^ e dz  ; 

Z OJ  . 

0 

Now  we  may  write  the  Fourier  transform  of  with  respect  to  0 as, 


^ e ^ d<p  = j . /3(z,y)  ^ (y)  cos  (i'0)l  e d0 


C|QINB  , BLaNK-NDT  fllMiD 
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wnich  as  the  series  is  uniformly  convergent  can  be  written. 


■I  3n 


^(z,0) 


= /3(z,y)  ^ o^^(y)  j cos  (v<p)  e 


■i  o).4> 
<P 


-i  u <t> 

(vtp)  e ^ d<t>  = n l&  (w  - r)  6 (w  + r){ 


hence 


^ j e d<i>  = /?(z,y)7r  (y)  + 5(o;^+i/;f 

A4 

This  then  is  the  Fourier  transform  with  respect  to  4>  of  which  we  may  now 

use  in  the  evaluation  of  equation  A3. 


Hence, 


+ o»  + «• 


T(p,y,<#>^)  = j T(z,w^)  ^(z,y)7T  ^ a^(y)  {5(w^-v)  + S{ui^*v)l 


xu>  A 

e dw,  dz 

<P 


+ o»  +0 


T(p,y,0^ ) 


^ j ^(z,y)  T(z,r)  dz  a^(y)  e^ 


This  then  gives  the  farfield  pressure  T(p,y,<^^)  in  terms  of  the  Fourier  trans- 
formed,  in  <t>,  nearfield  pressure  T(z,p). 

Sampling  of  the  Neafield  Pressure 

5.  Suppose  that  the  nearfield  pressure  T(z,<^)  is  not  known  as  a continuous 

function  of  z and  ({>  but  is  sampled  at  discrete  points  z,  ^ apart  then  we  may 
write  the  s ampled  nearfield  as 


t(z,<^)  = 


P=—  q=-“ 


T(pz,  q^)  8(z  - pz)  B{i>  - (^) 


and  hence  its  Fourier  transform  as,  (ref  5) 


•fo*  +0 


p=—  q=-. 
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Thus  we  may  rewite  equation  A5  as 


•f«  +«* 


I i /3(z»y)  dz  a^(y) 


which  by  expanding  ^(z>y)  and  letting  k = k cos  y may  be  written, 


T(p,y,<^^) 


^ • 


a /%  ' 


27T  ap 


and  as  x(k^,o)^)  may  be  written 


+ <»  +o 


^ I r ^ 


p=“-  q=“* 


we  have. 


't(pfy.'#>^)  = 


+ <«  +0.  +0 


2 — 
27t  ap0z 


ii  I I I f 


— p=**o* 


or  substituting  for  a^Cy)  and  letting  ka  sin  (y)  be  rewritten  as  (k^  - k 
we  have. 


2v  apz  4>  z ^ 4>  H^((kV^y/^a) 


l>=-~  p=-0.  q::_. 


Now  this^may  be  separated  into  the  correct  field  T(p,y,0^)  and  the  error  term 
Tg(p,y,0  ) as, 

^(p»yV)  = T(p,yV)  + Tg(p,yV)  Aj 

where,  the  correct  field  is  given  by, 

TCPjy,*^!  = — 2 ^ T(k  ,«.)  a (y)  ^ 

2ir^  ap  ^ ^ •' 


and  the  error  term  is  given  by, 


14. 


Tg(p,y,0^ ) 


^(ikp) 

2 

2it  ap  (f>z 


♦ o*  ( + 


II  I 

v=~^  [ p=-~  q=- 

(P  / o) 


a=—  2 


T(k,,».--2?)  «„(v) 

4> 


q=-~  " J 

(q/o) 

This  error  term  wi]l  be  negligible  as  long  as  the  following  inequalities  are 
obeyed, 

2tt 

— > k + k AlC 

- o zmax 
z 

27T  f 

where  k = and  k is  the  highest  axial  spatial  wave  number  present,  and 

C 7.tTif\.TC 


> M ka  sin  y + C + w 


where  is  the  highest  circumferential  spatial  wave  number  and  M and  C are 

constants  depending  on  the  accuracy  required,  which  for  most  practical  purposes 
take  ti'.e  values  M = 1 .4  and  C = ^6, 

Asymptotic  limits  of  the  circumferential  samplinc;  Criteria 


4.  It  is  interesting  to  consider  the  circumferential s ampling  criteria  for  the 
limit  of  large  radius  a.  Let  equations  AlO  and  A11  be  rewritten 


, I.  zmax  ^ 

k(i  ^ __) 


(Mkasiny+C)(^1  + yr ^ 


Now  let 


T^. — ^ = A a constant,  then 
[M  ka  sin  y * C)  ’ 


we  have  z < 


1 + A 


(1  + A)  sin  y (M  + 


Now  for  large  a,  ^ will  tend  to  zero  and  M tenda  to  1 (see  reference  b 

equation  9.2.3,  and  equation  A2)  and  the  limit  becomes 

* ^ (l  + d)  sin  y 


which  for  the  equitorial  plan  y = 90°  is  the  same  criterion  as  for  the  axial 


FIG.  3b  FOURIER  TRANSFORM  OF  THE  NORMAL  DERIVATIVE 
OF  THE  GREENS  FUNCTION  FOR  A CYLINDER 


i.Sb  NEAR  FIELD  PRESSURE  AND  IT  S FOURIER  TRANSFORM 


(y)  = 25.8 


FIS.  6a.  PREDICTED  FARPIBLD  FOR  VARIOUS  SAMPLING  CRITFEIA 


and 


FIG.  7b 


■©- 

3 


<K 


o 


U- 


dB 


-lo 


-30 


-30 


-40 


-SO 


■60 


I- 6 A 0'8X 


0 5X  = SAMPLE  SPACING 


O 


16 


32 


5 1 (V  ) 


i «<#, 


CIRCUMFERENTIAL  WAVENUMBER 


FIG.  7b  FOURIER  TRANSFORM  OF  SAMPLED 
NEARFIELD  (OMNIDIRECTIONAL) 


